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Abstract 

The behavior of the angular velocity of a test particle moving in a stable circular orbit in the 
CJ ' vacuum on the brane is considered. In the brane world scenario, the four dimensional effective 

Einstein equation acquire extra terms, called dark radiation and dark pressure, respectively, which 



X 



arise from the embedding of the 3-brane in the bulk. A large number of independent observations 
have shown that the rotational velocities of test particles gravitating around galaxies tend, as a 



> 

' function of the distance from the galactic center, toward constant values. By assuming a constant 

' tangential velocity, the general solution of the vacuum gravitational field equations on the brane 



can be obtained in an exact analytic form. This allows us to obtain the explicit form of the 



in 
o 

Q I projections of the bulk Weyl tensor on the brane, and the equation of state of the dark pressure 



as a function of the dark radiation. The physical and geometrical quantities are expressed in 
terms of observable/measurable parameters, like the tangential velocity, the baryonic mass and the 
radius of the galaxy. We also analyze the dynamics of test particles by using methods from the 
qualitative analysis of dynamical systems, by assuming a simple linear equation of state for the 
dark pressure. The obtained results provide a theoretical framework for the observational testing 
at the extra-galactic scale of the predictions of the brane world models. 
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I. INTRODUCTION 



The problem of the dark matter is a long standing problem in modern astrophysics. Two 
important observational issues, the behavior of the galactic rotation curves and the mass 
discrepancy in clusters of galaxies led to the necessity of considering the existence of dark 
matter at a galactic and extra-galactic scale. 

n 

The rotation curves of spiral galaxies jil!] are one of the best evidences showing the 
problems Newtonian gravity and/or standard general relativity have to face on the galac- 
tic/intergalactic scale. In these galaxies neutral hydrogen clouds are observed at large dis- 
tances from the center, much beyond the extent of the luminous matter. Since these clouds 
move in circular orbits with velocity Vtg{r), the orbits are maintained by the balance be- 
tween the centrifugal acceleration v^g/r and the gravitational attraction GM{r)/r'^ of the 
total mass M(r) contained within the orbit. This allows the expression of the mass profile 
of the galaxy in the form M(r) = rv^g/G. 

Observations show that the rotational velocities increase near the center of the galaxy, in 
agreement with the theory, but then remain nearly constant at a value of Vtnoo ~ 200 — 300 

n 

km/s [1], which leads to a mass profile M(r) = rvfg^/G. Consequently, the mass within a 
distance r from the center of the galaxy increases linearly with r, even at large distances 
where very little luminous matter has been detected. 

The second astrophysical evidence for dark matter comes from the study of the clusters 
of galaxies. The total mass of a cluster can be estimated in two ways. Knowing the motions 
of its member galaxies, the virial theorem gives one estimate, Myr, say. The second is 
obtained by separately estimating the mass of each individual members, and summing these 
masses, to give a total baryonic mass Mb- Almost without exception it is found that Myr is 
considerably greater than Mb, Mvt > Mb, typical values of Mvt/Mb being about 20-30 

This behavior of the galactic rotation curves and of the virial mass of galaxy clusters is 
usually explained by postulating the existence of some dark (invisible) matter, distributed 
in a spherical halo around the galaxies. The dark matter is assumed to be a cold, pressure- 
less medium. There are many possible candidates for dark matter, the most popular ones 
being the weakly interacting massive particles (WIMP) (for a review of the particle physics 
aspects of the dark matter see 2!]). Their interaction cross sections with normal baryonic 
matter, while extremely small, are expected to be non-zero, and we may expect to detect 
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them directly. 

It has also been suggested that the dark matter in the Universe might be composed of 
super-heavy particles, with mass > 10^° GeV. But observational results show that the dark 
matter can be composed of super-heavy particles only if these interact weakly with normal 
matter, or if their mass is above 10^^ GeV jsl- Scalar fields or other long range coherent 
fields coupled to gravity have also intensively been used to model galactic dark matter {4]. 

A general analysis of the possibility of an alternative four-dimensional gravity theory ex- 
plaining the dynamics of galactic systems without dark matter was performed by Zhytnikov 
and Nester jsj. From very general assumptions about the structure of a relativistic gravity 
theory (the theory is metric, and invariant under general coordinates transformation, has a 
good linear approximation, it does not possess any unusual gauge freedom and it is not a 
higher derivative gravity) a general expression for the metric to order [vjcf' has been de- 
rived. This allows to compare the predictions of the theory with various experimental data: 
the Newtonian limit, light deflection and retardation, rotation of galaxies and gravitational 
lensing. The general conclusion of this study is that the possibility for any four- dimensional 
gravity theory to explain the behavior of galaxies without dark matter is rather improbable. 

However, up to now no non-gravitational evidence for dark matter has been found, and 
no direct evidence or annihilation radiation from it has been observed yet. 

Therefore, it seems that the possibility that Einstein's (and the Newtonian) theory of 
gravity breaks down at the scale of galaxies cannot be excluded a 'priori. Several theoretical 
models, based on a modification of Newton's law or of general relativity, have been proposed 
to explain the behavior of the galactic rotation curves 6|. 

The idea of embedding our Universe in a higher dimensional space has attracted a con- 

n 

siderable interest recently, due to the proposal by Randall and Sundrum \ T\ that our four- 
dimensional (4d) spacetime is a three-brane, embedded in a 5d spacetime (the bulk). Ac- 
cording to the brane world scenario, the physical fields (electromagnetic, Yang-Mills etc.) 
in our 4d Universe are confined to the three brane. Only gravity can freely propagate in 
both the brane and bulk spacetimes, with the gravitational self-couplings not significantly 
modified. Even if the fifth dimension is uncompactified, standard 4d gravity is reproduced 
on the brane. Hence this model allows the presence of large, or even infinite non-compact 
extra dimensions. Our brane is identified to a domain wall in a 5d anti-de Sitter spacetime. 
For a review of the dynamics and geometry of brane Universes, see e.g. jsl. 
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Due to the correction terms coming from the extra dimensions, significant deviations 
from the standard Einstein theory occur in brane world models at very high energies j^. 
Gravity is largely modified at the electro-weak scale of about 1 TeV. The cosmological and 
astrophysical implications of the brane world theories have been extensively investigated in 
the physical literature Q, . 

The static vacuum gravitational field equations on the brane depend on the generally 
unknown Weyl stresses, which can be expressed in terms of two functions, called the dark 
radiation U and the dark pressure P terms (the projections of the Weyl curvature of the 



bulk, generating non-local brane stresses) |8|, Il3l |. 



Several classes of spherically symmetric solutions of the static gravitational field equations 



in the vacuum on the brane have been obtained in 
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18|. As a possible physical 



application of these solutions the behavior of the angular velocity Vtg of the test particles in 
stable circular orbits has been considered [l7,[l^. The general form of the solution, together 
with two constants of integration, uniquely determines the rotational velocity of the particle. 
In the limit of large radial distances, and for a particular set of values of the integration 
constants the angular velocity tends to a constant value. This behavior is typical for massive 
particles (hydrogen clouds) outside galaxies and is usually explained by postulating the 
existence of the dark matter. 

Thus, the rotational galactic curves can be naturally explained in brane world models, 
without introducing any additional hypothesis. The galaxy is embedded in a modified, 
spherically symmetric geometry, generated by the non-zero contribution of the Weyl tensor 
from the bulk. The extra terms, which can be described in terms of the dark radiation term 
U and the dark pressure term P, act as a "matter" distribution outside the galaxy. The 
particles moving in this geometry feel the gravitational effects of U, which can be expressed 
in terms of an equivalent mass (the dark mass) . The dark mass is linearly increasing with 
the distance, and proportional to the baryonic mass of the galaxy, Mu{r) ^ MB{r/ro) 17|. 

The exact galactic metric, the dark radiation and the dark pres sure in the flat rotation 
curves region in the brane world scenario has been obtained in 18|. 

Similar interpretations of the dark matter as bulk effects have been also considered in [Ji 
where it was also shown that it is possible to model the X-ray proflles of clusters of galaxies, 
without the need for dark matter. 

In the brane world scenario, the fundamental scale of gravity is not the Planck scale, but 
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another scale which may be at the TeV level. The gravitons propagating through the bulk 
space give rise to a Kaluza-Klein tower of massive gravitons on the brane. These gravitons 
couple to the energy-momentum term of the standard model fields, and could be produced 
under the appropriate circumstances as real or virtual particles. 

Another important effect that is expected in the brane world models is the presence 
of brane fluctuations, since rigid objects do not exist in the relativistic theory. It was 
proposed that, in the context of brane world scenarios with low tension r = Z"^, massive 

20 . The different 



brane fluctuations are natural dark matter candidates, called branons 
possibilities for branons as dark matter candidates, the parameter region in which branons 
behave as collisionless thermal relics (WIMPs), either cold or hot (warm) together with less 
standard scenarios in which they are strongly self interacting or produced non-thermally, 
and the possibilities of obtaining branons in hadron accelerators have been considered in 



2l| . Another possibility to test these models is through the annihilation of branons into 



photon or electron-positron pairs in the galactic halos, which could possibly be detected by 
gamma-ray or anti-matter detectors. 

A possible brane world model, which tries to give a justification for the hypothesis of the 
scalar field origin for the dark matter was proposed by 22| . The model contains two branes, 
on one of the branes lives the matter of the standard model of particles but on the other 
one, only spin-0 fundamental interactions are present. In the model the spin-0 fields are the 
dark matter. The scalar field contains an effective pressure, which avoids the collapse of the 
field fluctuations, implying that scalar field dominated objects, like galaxies, may contain a 
flat density profile in the center. 

Gravitational lensing and the study of the light deflection by black holes and galaxies 
is an important physical effect that could provide specific signatures for testing the brane 
world models (for a review of the gravitational lensing by brane world black holes see 23|). 
Observables related to the relativistic images of strong field gravitational lensing could in 
principle be used to distinguish between different brane world black hole metrics in future 
observations. The strong field limit approach was used in 2^] to investigate the gravita- 
tional lensing properties of brane world black holes, and the lensing observables for some 
candidate brane world black hole metrics have been compared with those for the standard 
Schwarzschild case. 

Brane World black holes could have significantly different lensing observational signatures 
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as compared to the Schwarzschild black holes. The bending angle predicted by the brane 
world models is much larger than that predicted by standard general relativistic and dark 
241 ] ■ The deflection of photons and the bendin g an gle of light in the 



matter models 



23 



constant tangential velocity region on the brane was considered in |18l . The bending 
angle predicted by the brane world models is much larger than that predicted by standard 
general relativistic and dark matter models. Therefore, the study of the gravitational lensing 
could discriminate between the different dynamical laws proposed to model the motion of 
particles at the galactic level and the standard dark matter models. 

Generally, the vacuum fleld equations on the brane can be reduced to a system of two 
ordinary differential equations, which describe all the geometric properties of the vacuum as 
functions of the dark pressure and dark radiation terms [16] . In order to close the system 
of vacuum fleld equations on the brane a functional relation between these two quantities is 
necessary. 

Hence, a flrst possible approach to the study of the vacuum brane consists in adopting 
an explicit equation of state for the dark pressure as a function of the dark radiation. A 
second method consist in closing the fleld equations on the brane by imposing the condition 
of the constancy of the rotational velocity curves for particle in stable orbits. 

It is the purpose of the present paper to consider the general behavior of the vacuum 
gravitational fleld equations in the brane world model in the region of constant tangential 
rotational velocity of test particles in stable circular orbits. Physically, this situation is 
characteristic for particles gravitating in circular orbits around the galactic center As 
a flrst step in our study we derive, under the assumption of spherical symmetry, the basic 
equations describing the structure of the vacuum on the brane, and the equation giving the 
behavior of the tangential velocity of the test particles as functions of the dark radiation 
and of the dark pressure. 

To obtain the tangential velocity of test particles we analyze the motion in an effective 
potential, also containing the effects of the bulk. By assuming that the brane is a flxed 
point of the bulk, we derive the equation giving the tangential velocity as a function of the 
00 component of the metric tensor only. In the region of constant tangential velocities, the 
general solutions of the gravitational equations can be obtained in an exact analytic form. 

This allows us to obtain the explicit form of the projections of the bulk Weyl tensor on 
the brane (the dark radiation and the dark pressure), as well as the components of the metric 
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tensor. The mass associated to the dark radiation (the dark mass) has a similar behavior as 
the dark matter at the galactic scale - it is a linearly increasing function of the distance, and 
it is proportional to the square of the tangential velocity. This result strongly supports the 
interpretation of the usual "dark matter" as a bulk effect. All the physical and geometrical 
quantities in our model are expressed in terms of observable/measurable parameters, like 
the tangential velocity, the baryonic mass and the radius of the galaxy. Some astrophysical 
applications of our model, including the possibilities of its observational testing, are also 
briefly considered. 

Since generally the structure equations of the vacuum on the brane cannot be solved 
exactly, we shall analyze the dynamics of test particles by using methods from the qualitative 
analysis of dynamical systems. In particular, we will investigate the general behavior of the 
tangential velocity on the brane by assuming a simple linear equation of state for the dark 
pressure. In this case generally we cannot reproduce the observed behavior of the test 
particles. 

The present paper is organized as follows. The field equations for the vacuum and the 
tangential velocity of a test particle on the brane are written down in Sections HTl and Unl The 
general solution of the field equations in the constant tangential velocity region is obtained 
in Section IIV[ The qualitative analysis of the field equations for a linear equation of state 
for the dark pressure is performed in Section |Vl We discuss and conclude our results in 
Section IVll 

II. THE FIELD EQUATIONS AND THE TANGENTIAL VELOCITY OF TEST 
PARTICLES FOR STATIC, SPHERICALLY SYMMETRIC VACUUM BRANES 

A. The field equations in the brane world models 

We start by considering a five dimensional (5d) spacetime (the bulk), with a single four- 
dimensional (4d) brane, on which matter is confined, only gravity can probe the extra 
dimensions. The 4d brane world {^^^M,g^,^) is located at a hypersurface [B {X^) = O) in 
the 5d bulk spacetime {^^^M, gAs), of which coordinates are described by X"^, A = 0,1, ... ,4. 
The induced 4d coordinates on the brane are x'^, ii = 0,1, 2, 3. 
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The action of the system is given by 



•S" — Sbuik + Si 



''brane i 



(1) 



where 



and 



hulk 



(5)M 



2kl 



(FX. 



S, 



brane 



(4)M 



-'-'brane 

L^5 



d X, 



(2) 



(3) 



where = SttG^ is the 5d gravitational constant, ^^^R and ^^^Lm are the 5d scalar curvature 
and the matter Lagrangian in the bulk, L^rane {dais, i^) is the 4d Lagrangian, which is given 
by a generic functional of the brane metric gap and of the matter fields ip, is the trace 
of the extrinsic curvature on either side of the brane, and A5 and (the constant brane 
tension) are the negative vacuum energy densities in the bulk and on the brane, respectively. 
The Einstein field equations in the bulk are given by 



(4) 



where 



(5)t . - o ^''^^^"^ I (5) (5)r (r\ 

"~ 5(5) glJ^ ^' 

is the energy-momentum tensor of bulk matter fields, while T^j, is the energy-momentum 
tensor localized on the brane and which is defined by 



T, 



6Ll)rane j 

~ ~r gpLV ^brane- 



Sg 



(6) 



The delta function 6 {B) denotes the localization of brane contribution. In the 5d space- 
time a brane is a fixed point of the Z2 symmetry. The basic equations on the brane are 
obtained by projections onto the brane world. The induced 4d metric is gjj = ^^^gij — njUj, 
where nj is the space-like unit vector field normal to the brane hypersurface ^^^M. In the 
following we assume ^^'^Lm = 0. 

Assuming a metric of the form ds"^ = {niUj + gij)dx^dx'^ , with riidx^ = dx the unit 
normal to the x = constant hypersurfaces and gu the induced metric on x = constant 
hypersurfaces, the effective 4d gravitational equation on the brane (the Gauss equation), 

9I: 



takes the form 



G 



111/ 



AQuu + k^T^jy -\- k^S^y 



E 



(7) 



where S^i, is the local quadratic energy-momentum correction 

11 1 

and Efj^^ is the non-local effect from the free bulk gravitational field, the transmitted pro- 
jection of the bulk Weyl tensor Cjajb, Ejj = CjAjsn'^n^ , with the property Eu — >■ 
Efiu^i^j as X 0- We have also denoted fc| = SttG, with G the usual 4d gravitational 
constant. 

The 4d cosmological constant, A, and the 4d coupling constant, k^, are related by A = 
^5(^5 + k'^^l/Q)/2 and = k^Xb/6, respectively. In the limit A^^ — we recover standard 
general relativity 9j. 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation 
of the energy- momentum tensor of the matter on the brane, D^T^'^ = 0, where D,y denotes 
the brane covariant derivative. Moreover, from the contracted Bianchi identities on the 
brane it follows that the projected Weyl tensor obeys the constraint D^E^" = k^DyS^" . 

The symmetry properties of E^j^ imply that in general we can decompose it irreducibly 
with respect to a chosen 4-velocity field u'^ as [l^ 



Efxu — —k 



(9) 



where k = k^/k^, h^y = g^^ + u^Uy projects orthogonal to u^, the "dark radiation" 
term U = —k'^E^^u^u^ is a scalar, = k~^h^Eaj3U^ is a spatial vector and P^y = 
—k~* \h(^^ "hy^ ^ — Ea(3 is a spatial, symmetric and trace-free tensor. 

In the case of the vacuum state we have p = p = 0, T^^ = 0, and consequently 5"^,^ = 0. 
Therefore the field equation describing a static brane takes the form 

Rf,^ = -E^^ + A^^^, (10) 

with the trace R of the Ricci tensor R^^ satisfying the condition R = R^ = 4A. 

In the vacuum case E^^, satisfies the constraint D^E^^ = 0. In an inertial frame at any 
point on the brane we have = 6q and huu = diag(0, 1, 1, 1). In a static vacuum = 
and the constraint for E^^i, takes the form [8j 

^D^U + ^UA^ + D-'P^, + A'^P,, = 0, (11) 

where A^ = u^Di^u^ is the 4-acceleration. In the static spherically symmetric case we 
may chose A^ = y4(r)r^ and P^,y = P(r) (r^r^, — \h^y), where A{r) and P{r) (the "dark 



pressure" although the name dark anisotropic stress might be more appropriate) are some 
scalar functions of the radial distance r, and is a unit radial vector 13|. 



B. The motion of particles in stable circular orbits on the brane 

In brane world models test particles are confined to the brane. Mathematically, this 
means that the equations governing the motion are the standard 4d geodesic equations 
HQ . However, the bulk has an effect on the motion of the test particles on the brane 
via the metric. Since the projected Weyl tensor effectively serves as an additional matter 
source, the metric is affected by these bulk effects, and so are the geodesic equations. This 
is in contrast to Kaluza-Klein theories, where matter travels on 5d geodesies. 

In order to obtain results which are relevant to the galactic dynamics, in the following 
we will restrict our study to the static and spherically symmetric metric given by 

ds^ = -e'^'Ut" + e'^^'Ur^ + r^c/fi^, (12) 

where dVL^ = dO'^ + sin^ OdcfP'. 

The Lagrangian C for a massive test particle traveling on the brane reads 

£=l(^-e-i2^eV + r2fi2), (13) 

where the dot means differentiation with respect to the afl&ne parameter. 

Since the metric tensor coefficients do not explicitly depend on t and f2, the La- 
grangian f|T3|) yields the following conserved quantities (generalized momenta): 

- e^'^'H = E, r^Q = L, (14) 

where E is related to the total energy of the particle and L to the total angular momentum. 
With the use of conserved quantities we obtain from Eq. (|T3|) the geodesic equation for 
massive particles in the form 

e-V^ + (l + ^) = E\ (15) 

The second term of the left-hand side can, in some cases, be interpreted as an effective 
potential. For instance, for the Schwarzschild space-time, where e'^'^^ = 1, the kinetic term 
is position independent. In that case the notion of an effective potential is appropriate. In 
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other cases, even one can still compute the turning points of the kinetic term, however, the 
effective potential interpretation is lost. 

For the case of the motion of particles in circular and stable orbits the 'potential' must 
satisfy the following conditions: a) r = (circular motion) b) dVf.ff/dr = (extreme motion) 
and c) d'^Veff/dr ^ \^^^^ > (stable orbit), respectively. Conditions a) and b) immediately 
give the conserved quantities as 

E' = e^(l + ^). (16) 



and 



respectively. Equivalently, these two equations can be rewritten as 

= = (18) 

l-rv'/2 l-rv'/2 ^ ^ 

We define the tangential velocity vtg of a test particle on the brane, as measured in terms 
of the proper time, that is, by an observer located at the given point, as 26| 

By using the constants of motion, and eliminating the quantity L, we obtain the expres- 
sion of the tangential velocity of a test particle in a stable circular orbit 251] on the brane 

as 

< = -Y- (20) 

Let us emphasize again that the function v' is obtained by solving the field equations 
containing the bulk effects as additional matter terms; we consider this in Section IIIII 



C. The gravitational field equations for a static spherically symmetric brane 



With the metric given by ( |T2l) the gravitational field equations and the effective energy- 
momentum tensor conservation equation in the vacuum take the form jl6. 17] 



— e 



-A 



1 A' 



1 



z/' 1 \ 1 



r2 k^X^ 



f/ + A, 



(f/ + 2P) - A, 



(21) 
(22) 
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2f/ + P r(2f/ + P)' 
where ' = djdr. In the following we denote 

167rG . . 

As for the motion of the test particle on the brane we assume that they follow stable 
circular orbits, with tangential velocities given by Eq. fl20|) . Thus, the rotational velocity of 
the test body is determined by the metric coefficient exp [y) only. 

The field equations ( l2Ti) -( l22l) yield the following effective energy density, radial and or- 
thogonal pressure 

p^fl^ = 3af/ + A, (26) 
P^ff = af/ + 2aP-A, (27) 
Pf = aU-aP-A, (28) 

which by rewriting become 

pcS _ pes _ 2pcfT = 4^ = const. (29) 

This is expected for the 'radiation' like source, the projection of the bulk Weyl tensor, 
which is trace-less Ei^ = 0. 

III. STRUCTURE EQUATIONS OF THE VACUUM IN THE BRANE WORLD 
MODELS 

Eq. (12T1) can immediately be integrated to give 

e- = 1 - ^1 - ^M-M - (30) 

r r 3 

where Ci is an arbitrary constant of integration, and we denoted 

pr 

GMu (r) = 3a / U{rydr. (31) 
Jo 

The function Mu is the gravitational mass corresponding to the dark radiation term 
(the dark mass). For U = the metric coefficient given by Eq. (I30l) must tend to the 
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standard general relativistic Schwarzschild metric coefficient, which gives Ci = 2GM, where 
M = constant is the baryonic (usual) mass of the gravitating system. 

By substituting u' given by Eq. into Eq. (12^ and with the use of Eq. (I5U]) we obtain 
the following system of differential equations satisfied by the dark radiation term U, the 
dark pressure P and the dark mass M[/, describing the vacuum gravitational field, exterior 
to a massive body, in the brane world model 16l |: 

dU _ {2U + P) [2GM + GMu + a{U + 2P) r^] - fAr^ dP 6P 

^2 (1 _ 2GM _ _ A^2) ^ rf7 ~ V' ^^^^ 

V r r 3 / 

dMu 3a 2 



dr a 

The tangential velocity can be written as 



-r'U. (33) 



2 1 2GM + GMu + a{U + 2P) - fAr^ 

V- — 



(34) 



2 r (l - ^"^^^ - "^^-^c/ _ A^2) 

V r r 3 / 

To close the system a supplementary functional relation between one of the unknowns U, 
P, Mjj and Vtg is needed. Once this relation is known, Eqs. (I5^ - (I2U]) give a full description 
of the geometrical properties and of the motion of the particles on the brane. 

The system of equations fl32|) and fl33|) can be transformed to an autonomous system of 
differential equations by means of the transformations 

2GM GMu 2 o 2,,^o 2A 

q = \ 1 r , /i = 3ar U + 3r A, 

r r 3 

p = Sar^P - 3r^A, 9 = \nr. (35) 

We shall call fi and p the "reduced" dark radiation and pressure, respectively. 

With the use of the new variables given by Eqs. (135|1 . Eqs. ( l32i) and ( 133|) become 

^^_(2,+ri[,+ i(, + 2p)]_^^^ (37) 

Eqs. ( 132!) and (133|) . or, equivalently, (!36|) and (1371) may be called the structure equations 
of the vacuum on the brane. In order to close this system an "equation of state", relating 
the reduced dark radiation and the dark pressure terms is needed. Generally, this equation 
of state is given in the form P = P{U). 
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In the new variables the tangential velocity of a particle in a stable circular orbit on the 
brane is given by 

2 _ 1 g + |(/^ + 2p) 

""^^ " 2 1-q ■ 
By using the expression of the tangential velocity, Eq. fl371) can be rewritten as 



(38) 



^ (/i + 2p) = -2 (2/i + p) vl + 2/i - 2p. 



(39) 



Eq. (1391) allows the easy check of the physical consistency of some simple equations of 
state for the dark pressure. The equation of state /i + 2p = immediately gives v'fg = 1, 
implying that all test particle in stable circular motion on the brane move at the speed of 
light, a fact that is contradicted by the observations at the galactic scale. The equation of 
state 2/i + p = gives /i = /Uq/t^, where /iq = constant is an arbitrary integration constant, 
U = /io/3ar^ and GMu = —fio/3r^, respectively. In the limit of large r, the tangential 
velocity vfg tends to zero, Vtg 0. Therefore, this model seems also to be ruled out by 
observations. The case fi = p gives fi{9) = fioexp [—2 J v^g{9)d9~\, = constant, and 
q{e) = {2vl-^^)/{l + 2vl). 

The dark radiation and the dark pressure can be obtained as a functions of the tangential 
velocity in a closed analytical form in the important case of a linear equation of state of the 
form 

p=(r-2)/i + A (40) 
with r and (3 arbitrary constants. Then the reduced dark radiation can be obtained as 

2r 



/,(^^) = ^^2(3-r)/(2r-3)g^p 
X <^ Co 



2r-3 



2r-3 



< id) d9 

,-2(3-r)/(2r-3) 



exp 



2r 



2r - 3 



vl (6) de 



(41) 



where Cq is an arbitrary integration constant. Hence, if the velocity profile of a test particle 
in stable circular motion is known, one can obtain all the relevant physical parameters for a 
static spherically symmetric system on the brane. 



IV. DARK RADIATION, DARK PRESSURE AND THE GALACTIC METRIC 
IN THE CONSTANT TANGENTIAL VELOCITY REGION 

The galactic rotation curves provide the most direct method of analyzing the gravitational 
field inside a spiral galaxy. The rotation curves have been determined for a great number of 
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spiral galaxies. They are obtained by measuring the frequency shifts z of the light emitted 
from stars and from the 21cm radiation emission from the neutral gas clouds. Usually, the 
astronomers report the resulting z in terms of a velocity field Vtg. The observations show 
that at distances large enough from the galactic center 

vtg ^ 200 - 300km/s = constant. (42) 

This behavior has been observed for a large number of galaxies 

In the following we use this observational constraint to reconstruct the metric of a galaxy 
on the brane. The constancy of Vtg allows us to express the function q from Eq. (l20l) as 



giving immediately 

dQ ^ 3(l + 2t;y V^^^ ' "c^^ 

From Eq. fl39l) we obtain 



d(l 1 (d^ , 2^') (44) 



* + 2|^-2(l-2,g,-2(l + „?>, (45) 
which can be used to simplify Eq. ( 14^ into 

By eliminating dq/dO between Eqs. ( l36l) . and ( 1371) and by using the expression of q given 
by Eq. fj^3l) . we obtain the equation of state for the dark pressure in the constant tangential 
velocity region as 

p = j^ii - 3. (47) 

This shows that, as a function of the "reduced" dark radiation, in the coordinate 9, the 
"reduced" dark pressure obeys a generalized linear equation of state. 

By using the generalized linear equation of state, Eq. (1371) can be written as 

dii 

— = -n/i + m, (48) 



where 

2vl \^ + v^ 



n — 



3 (2 + 



+ 3vl + 2 



„,2 ' '^^tg 
Vtg 



(49) 
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and 



m 



2 

2 + < 



Eq. fl48|) can be immediately integrated to give 



m 



n 



(50) 



(51) 



with C an arbitrary constant of integration. 

In the standard radial coordinate r the "reduced" dark radiation energy density is given 

by 



n(r) = — + ^~ ' "^g^ *g 



while the "reduced" dark pressure is 

p{r) - 



3 + < C 

Vtg r 



(52) 



(53) 



From these equations we find the variation of the dark radiation and dark pressure in the 
constant tangential velocity region as 



and 



U{r) = — 



+ 



3 + < C 
Vtg r 



2u 



tg 



1 + vig r2 



-3A 



+ 3A 



(54) 



(55) 



respectively. 

The metric tensor component exp(z/) can be found from the expression of the tangential 
velocity as 

(56) 



with 6*2 an arbitrary constant of integration, while exp(A) is given by 



2GM 



C 



'tg 



A 



l+Vl+Vfg 3 



(57) 



r {1 — n) r"^ 

The function M(/(r), which gives the mass associated with the dark radiation, is given 



by 



GMuir 



C 



(1 - n) r"-i ' 1 + v'f + vf 



+ 



r. 



(58) 



Hg -r '^tg 

Hence, the problem of the determination of the metric and of the dark radiation and dark 
pressure in the constant tangential velocity region around galaxies in the brane world model 
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is completely solved. It remains to determine the values of the integrations constants C and 
C2. 

To do this we assume that in the presence of matter the bulk effects on the brane are 
negligible small. This definitely represents a very good approximation for the galaxy, which 
has as only constituent "normal" matter, with a density much higher than the density of the 
matter outside the galaxy. On the other hand, at large distances from the galactic center 
the cosmological background dominates the dynamics. 

Therefore, by assuming that the radius of the baryonic matter distribution in the galaxy 
is we can match our solution iov r — R with the Schwarzschild-de Sitter metric, so that 

(59) 



„^ I 

e \r=R 



2GM A „2 

1 R^, 

R 3 ' 



and 



2GM A 



-R 



respectively. This gives 



and 



Co^R 



-2v: 



tg 



R 3 
2GM\ 



1 - 



(60) 
(61) 
(62) 

respectively. Therefore the dark mass, the dark radiation and the dark pressure terms can 



R 



C = -(l-n) 



tgj "tg 



1 



finally be written as 



aU{r) 



GMu{r) 



(1 + 



tg) "^"tg , 



1 + < + vt 



1 - 



'tg 



3 (1 + < + <) 



aP{r) = -1 



n] 



ig) 



(63) 
(64) 

(65) 



3 (1 + vl + 

Hence, the components of the Weyl tensor from the bulk can be obtained in terms of the 
observable quantities, like the baryonic radius of the galaxy and the tangential velocity of 
test particles gravitating in stable circular orbits around the galactic center. 

The brane world metric in the constant tangential velocity region has a singularity at some 
r = rh, which can be obtained by solving the non-linear algebraic equation exp(— A) = 0. 
Depending on the numerical values of n this equation may have one or more positive roots, 
which define an event horizon. Therefore the singularity at r = is hidden and cannot be 
seen by an external observer. 
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V. QUALITATIVE ANALYSIS OF THE STRUCTURE EQUATIONS OF THE 
VACUUM ON THE BRANE FOR GIVEN SIMPLE EQUATIONS OF STATE OF 
THE DARK PRESSURE 



In the present Section we investigate the general behavior of the tangential velocity on the 
brane by assuming a simple linear equation of state for the dark pressure. Since generally the 
structure equations of the vacuum on the brane cannot be solved exactly, we shall analyze 
the dynamics of test particles by using methods from the qualitative analysis of dynamical 
systems, by closely following the approach of 

We consider the case in which the dark pressure is proportional to the dark radiation 
P = 7f/, where 7 is an arbitrary constant, which can take both positive and negative values. 
In the reduced variables fi and p the linear equation of state is p = 7/2, and the structure 
equations of the gravitational field on the brane have the form 

fg=l^-1. (66) 

The tangential velocity of the particles in stable circular orbit is given by 



Let us firstly analyze the special case where 7 = —1/2. Then, by virtue if the second 
Eq. fl67j) . we obtain two possible solution, either g = yU = 0org = /i = 2/3. Since both 
solution are physically uninteresting, implying either Vtg = or Vtg = 1, let us assume 
henceforth that 7 7^ —1/2, and rewrite the system of equation into the following form 

^_ 2(1-7) 7 + 2 

d9 1 + 27^ 1 + 27 1-q ' ^'^ 

which is finally written as 

| = ^f + «. (71) 

where we have denoted 

^= I M , A=\ . _ , h ^=1 ...„l..l^„^ I (72) 





18 



The system of equations (17T1) has two critical points Xq = (0, 0) and = (3(1 — 
7)/(7^ + 7 + '''))3(1 — 7)/(7^ + 7 + 7)). For 7 = 1, the two critical points of the system 
coincide. Depending of the values of 7, these points lie in different regions of the phase space 
plane (g, /i). 

Since the term ||^||/||^|| as 1 1^| | ^ 0, the system of equations (|7Ti) can be linearized 
at the critical point Xq. The two eigenvalues of the matrix A are given by ri = —1 and 
r2 = 2(1 — 7)/(l + 27), and determine the characteristics of the critical point Xq. For 
7 G (—00,— 1/2) U [1,00) both eigenvalues are negative and unequal. Therefore, for such 
values of 7 the point Xq is an improper asymptotically stable node. 

If 7 G (—1/2,1), we find one positive and one negative eigenvalue, which corresponds 
to an unstable saddle point at the point Xq. Moreover, since the matrix dA/d^{XQ) has 
real non-vanishing eigenvalues, the point Xq is hyperbolic. This implies that the properties 
of the linearized system are also valid for the full non-linear system near the point Xq. It 
should be mentioned however, that this first critical point is the less interesting one from a 
physical point of view, since it corresponds to the 'trivial' case where both physical variables 
vanish. 

Before performing a similar analysis for the second critical point, it is worth mentioning 
that the structure equations can be solved exactly for the value 7 = —2. In that case, the 
non-linear term n in Eq. flTTj) identically vanishes, and the system of equations becomes a 
simple linear system of differential equations. For 7 = —2 the two eigenvalues of A are given 
by Ti = —1 and r2 = —2, the two corresponding eigenvectors are linearly independent and 
the general solution can be written as follows 



where go = Q'(O) and ^0 = Ai(0). One can easily transform this solution back into the radial 
coordinate r form by using 6 = log(r) and get 



This solution can be used to analyze the resulting tangential velocity profile as a function 
of the radius r, which turn out to have the same form as predicted by the Newtonian analysis. 
Therefore such an equation of state, namely p = — 2/i, cannot account for fiat rotation curves 
in galaxies, see the end of Section IIIII 




(73) 




(74) 
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Let us now analyze the qualitative behavior of the second critical point Xj. To do this, 
one has to Taylor expand the right-hand sides of Eqs. fl69l) and flTOj) around and obtain the 
matrix A which corresponds to the system, linearized around X^. This linearization is again 
allowed since the resulting non-linear term, n say, also satisfies the condition | |/| I^CI I ~^ 

1 1^1 1 — ^ Xy. The resulting matrix reads 

^ = I . ^ ^1 (75) 



^(2+7)2(1+27) 2+7^ 



and its two eigenvalues are given by 



-3 - 67 ± + 873 + 13272 - 287 - 47 

472 + IO7 + 4 



— o — w y s- \/ iw y -r 07- -r ±0^7 — ^07 — ^1 



For —0.5 < 7 < 0.674865 the argument of the square root becomes negative and the 
eigenvalues complex. Moreover, the values 7 = — 1/2 and 7 = —2 have to be excluded, since 
the eigenvalues in Eq. (175!) are not defined in these cases. However, both cases have been 
treated separately above. 

If 7 G {—00, — 1/2) U (1, 00), then X^ corresponds to an unstable saddle point and for 7 G 
(0.674865, 1) it corresponds to an asymptotically stable improper node. More interesting is 
the parameter range 7 G (—1/2,0.674865) where the eigenvalues become complex, however, 
their real parts are negative definite. Hence, for those values we find an asymptotically stable 
spiral point at X^. Since this point is also an hyperbolic point, the described properties are 
also valid for the non-linear system near that point, see Fig. [H This qualitative study of the 
behavior of the structure equations of the vacuum on the brane shows that not all equations 
of state of the dark pressure would predict flat rotation curves. 

With Eq. ( H7I) we were able to explicitly construct the equation of state that yields 
rotation curves of constant velocity. Since these velocities are about 200 — 300km/s, we 
obtain as an equation of state 

p ^ (3 - 7) X lOV - 3, (77) 

which is of the general form p = afi — b. 

The above qualitative analysis can easily be repeated for the physically relevant equation 
of state. The first difference is that Xq = (0, 0) is no longer a critical point to the system, and 
the critical point X^ gets a contribution due to the parameter b. Moreover, X^ is replaced 
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FIG. 1: The left figure shows the phase space plot fi = fi{q) of the system (f69]l ~ (i70]) for 7 = 
0.15. The dot represents the critical point and one clearly sees that corresponds to an 
asymptotically stable spiral point. In the right figure 7 = 0.70, and there the critical point 
corresponds to an asymptotically stable improper node. 

by a pair of critical points, related to the two square roots of a quadratic equation that 
reduces to a linear equations as a 0. The critical point in the p > and fi > region 
in phase space corresponds to an asymptotically stable improper node, and the qualitative 
picture is that of the right plot of Fig. [H 



VI. DISCUSSIONS AND FINAL REMARKS 



The galactic rotation curves and the mass distributions in clusters of galaxies continue 
to pose a challenge to present day physics. One would like to have a better understanding 
of some of the intriguing phenomena associated with them, like their universality and the 
very good correlation between the amount of dark matter and the luminous matter in the 
galaxy. To explain these observations, the most commonly considered models are based on 
particle physics in the framework of Newtonian gravity, or of some extensions of general 
relativity [6|]. 

In the present paper we have considered, and further developed, an alternative view to 
the dark matter problem 17|, llSl], namely, the possibility that the galactic rotation curves 
and the mass discrepancy in clusters of galaxies can naturally be explained in models in 
which our Universe is a brane in a multi- dimensional spacetime. The extra terms in the 
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gravitational field equations on the brane induce a supplementary gravitational interaction, 
which can account for the observed behavior of the galactic rotation curves. As one can see 
from Eq. (1631) . in the limit of large r ^ R, and by taking into account that v^g ^ 1, we obtain 
GMu ^ fj^^r, a behavior which is perfectly consistent with the observational data [1], and 
is usually attributed to the existence of the dark matter. Since inside the galaxy Mjj ^ 0, 
it follows that the constant tangential velocity of the test particles is determined by the 
baryonic mass Mb and the radius R of the galaxy by the equation v^^ ~ GMb/R- This 
gives for the dark mass the scaling relation 



Mu 



(78) 



Mb R 

By using the simple observational fact of the constancy of the galactic rotation curves, the 
galactic metric and the corresponding components of the projected Weyl tensor E^i, (dark 
radiation and dark pressure) can be completely reconstructed, without any supplementary 
assumption. 

Weak lensing provides a unique probe of the gravitational potential on large scales. Hence, 
in the context of dark matter, it can provide constrains on the extent and shapes of dark 
matter halos. Furthermore, it can test alternative theories of gravity (without dark matter). 
Due to the fixed form of the galactic metric on the brane, in our model the light bending 
angle is a function of the tangential velocity of particles in stable circular orbit and of the 
baryonic mass and radius of the galaxy. The specific form of the bending angle is determined 
by the brane galactic metric, and this form of the metric is different as compared to the 



other dark matter models 



on g-ra nge sel: 



-interacting scalar fields, MOND, non-symmetric 
gravity etc.). As shown in 18|, ll9|, l23|, l2^, the gravitational light deflection angle is much 
larger than the value predicted by the standard general relativistic approach. Even when 
we compare our results with standard dark matter models, like the isothermal dark matter 
halo model, we still may find significant differences in the lensing effect. Therefore the study 
of the gravitational lensing may provide the evidence for the existence (or non-existence) of 
the bulk effects on the brane. On the other hand, since in this model there is only baryonic 
matter, all the physical properties at the galactic level are determined by the amount of the 
luminous matter and its distribution. 



To test alternative theories of gravity by using lensing one can use two approaches 28 1. 
A measurement of the radial dependence of the lensing signal gives the best accuracy. This 
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method requires the knowledge of the deflection law. The most direct test is the detection of 
the azimuthal variation of the lensing signal around the lenses, and the measurements does 
not require knowledge of the deflection law. Any anisotropy in the lensing signal caused by 
the galaxy decreases as and therefore is negligible in the bulk effects dominated region. 

The detection of an anisotropy in the lensing signal around galaxies may provide the best 
test for discriminating between dark matter and alternative gravity models, including the 
present model. The results of a study of weak lensing by galaxies based on 45.5 deg^ of 
Rc band imaging data from the Red-Sequence Cluster Survey were presented in j29|. A 
significant flattening was detected, which implies that the halos are well aligned with the 
light distribution, and an isotropic lensing signal is excluded with 99.5% confidence. 

From an observational point of view an important problem is to estimate an upper bound 
for the cutoff of the constancy of the tangential velocities. The idea is to consider the point 
at which the decaying density profile of the dark radiation associated to the galaxy becomes 
smaller than the average energy density of the Universe, and the cosmological expansion 
dominates the dynamics of the particles. Let the value of the coordinate radius at the point 
where the two densities are equal be Ru- Then at this point 3a U{Ru) = {8'KG)puniv, where 
Puniv is the mean energy density of the universe. 

An alternative estimation of Ru can be obtained from the observational requirement that 
at the cosmological level the energy density of the dark matter (which in our case is the dark 
mass associated to the dark radiation) represents a fraction Qdai ~ 0.25 of the total energy 
density of the flat universe, with density parameter Q = 1. Therefore the "dark matter" 
contribution inside a ball of radius Ru is given by AnQDAiRuPcrit/^, which gives 

1 1 



-vtg, (79) 



v/2rW Hz) Ho 

where h{z) is the Hubble constant normalized to its local value: h'^{z) = Qm (1 + z)^ + Q^, 
Qm is the total mass density parameter and is the dark energy density parameter. For a 
tangential velocity of the order of 200 km/s it follows that Ru ~ 2.8 Mpc. Therefore, we pre- 
dict that the flat rotation curves region should extend up to distances of a few megaparsecs 
from the galactic center. 

If we assume that the flat rotation curves extend indefinitely, the resulting spacetime is 
not asymptotically flat, but of de Sitter type. This is due to the presence of the cosmological 
constant A on the brane. Observationally, the galactic rotation curves remain flat to the 
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farthest distances that can be observed. Therefore, from both the observational and the 
theoretical point of view it is important to estimate the possible role of the cosmological 
constant on the extra galactic dynamic. By assuming that the cosmological constant has a 
numerical value of the order of A ^ 3 x 10~^^ cm~^ [30|, at a distance of r = 100 kpc, for a 
galaxy with mass M = 10^^ Mq, the quantities GM/r ^ 4.94 x 10"^ and Ar^ ^ 2.7 x 10"^ are 
roughly of the same order of magnitude. For clusters of galaxies with masses of the order of 
10^^ Mq and radii of the order of 2 Mpc, we have GM/r ^ 2.47 x lO"*' and Ar^ ^ 1.08 x 10"^ 
respectively. 

Hence the presence of a non-zero cosmological constant may play an important role 
in the dynamics of the test particles in stable circular orbits at the outer boundary of 
the galaxies or clusters of galaxies, since it may generate physical effects that are of the 
same order of magnitude as the extra-dimensional effects from the bulk. Therefore, for the 
correct estimation in regions far away from the galactic center of the dark radiation, of the 
dark pressure, of the tangential velocities of test particles and of other related effects (like, 
for example, the bending angle of the light) one must take into account the effect of the 
cosmological constant on the brane. 

In the present model all the relevant physical quantities, including the dark mass (playing 
the role of the dark matter on the brane), the dark radiation and the dark pressure describing 
the non-local effects due to the gravitational field of the bulk, are expressed in terms of 
observable parameters - the tangential velocity, the baryonic mass and the radius of the 
galaxy. 

There are three observationally testable predictions of the model: a) the constant rotation 
curve region may extend much farther than presently observed, up to a few megaparsecs, 
b) the weak lensing effect in the dark radiation dominated region is different from that 
predicted by the standard dark matter models, and other alternative gravitational theories 
and c) the numerical value of the maximum rotational velocity is mainly determined by the 
"normal" matter content of the galaxy. 

As for the main advantages of our model, on the one hand it is esthetically attractive 
- "matter" is a manifestation of geometry, it is relatively simple and, on the other hand, 
it is directly testable observationally. The best testing ground for brane world models is 
the extra-galactic space, where the effects of galactic/extra-galactic matter perturbations 
are minimal, and where we expect that the bulk effects appear in a "pure" form, and they 
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dominate the dynamics of test particles (gas clouds). Moreover, the extra-galactic space is 
accessible to direct optical or radio observation. In the solar system the bulk effects are very 
small, and difficult to observe. 

Therefore, these results open the possibility of testing the brane world models by using 
direct astronomical and astrophysical observations at the galactic or extra galactic scale. 
In this paper we have provided some basic theoretical tools necessary for the in depth 
comparison of the predictions of brane world models and of the observational results. 
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